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Abstract 

We propose a CUSUM type test for constant correlation that goes beyond a 
previously suggested correlation constancy test by considering Spearman's rho in 
arbitrary dimensions. Since the new test does not require the existence of any mo- 
ments, the applicability on usually heavy-tailed financial data is greatly improved. 
We calculate the asymptotic null distribution using an invariance principle for the 
sequential empirical copula process. The limit distribution is free of nuisance pa- 
rameters and critical values can be obtained without bootstrap techniques. We give 
a local power result and analyze the test's behavior in small samples. 

Keywords: Copula, Mixing, Multivariate sequential empirical process, Robustness, 
Structural break 
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1. Introduction 



Recently, Wied, Kramer and Dehling (2012) proposed a fluctuation test for constant 
correlation based on the Bravais-Pearson correlation coefficient. The test, which will be 
referred to as BPC test in the following, is useful e.g. in financial econometrics when 
a practitioner wants to analyze if correlations of asset returns change in time, see e.g. 



Longin and Solnik (1995) and Krishan et al. (2009) for the relevance of this question. It 



complements former approaches by e.g. Galeano and Pena (2007) and Aue et al. (2009). 
However, one major drawback of this test is the fact that the limit distribution is derived 



under the condition of finite fourth moments (similar to Aue et al. , 2009 ). This is a critical 



assumption because the existence of fourth moments in usually heavy-tailed financial 



returns is doubtful, see e.g. Grabchak and Samorodnitsky (2010), Kramer (2002) and 



Amaral et al. (2000). 



This paper presents a fluctuation test for constant correlation based on Spearman's rho 
and the sample version of it which imposes no conditions on the existence of moments. 
There are several advantages of Spearman's rho compared to the Bravais-Pearson correla- 
tion: In many situations, e.g. if the data is non-elliptical, the Bravais-Pearson correlation 
may not be an appropriate measure for dependence. It is confined to measuring linear 
dependence, while the rank-based dependence measure Spearman's rho quantifies mono- 
tone dependence. If the second moments do not exist, the Bravais-Pearson correlation is 
not even defined, while Spearman's rho does not require any moments. 
Spearman's rho is probably the most common rank-based dependence measure in eco- 



nomic and social sciences, see e.g. Gaifiler and Schmid (2010), who propose tests for 



equality of rank correlations, and the references herein. In addition, Spearman's rho 



often performs better in terms of robustness than the Bravais-Pearson correlation. Em 



brechts et al. (2002 ) discuss several other pitfalls and possible problems for a risk manager 



who simply applies the Bravais-Pearson correlation. 

Therefore it is natural in the context of testing for changes in the dependence structure of 
random vectors to extend the BPC test to a test for constant Spearman's rho. As expected 



from the theory of dependence measures, this test is applicable in more situations: It has 
a much better behavior in the presence of outliers and there are no conditions on the 
existence of moments. In addition, the test is applicable in arbitrary dimensions, while 
the BPC test is designed for bivariate random vectors. Similarly to the BPC test, the test 



bases on successively calculated empirical correlation coefficients in the style of Ploberger 



et al. (1989), Lee et al. (2003) or Galeano (2007) and the limit distribution of our test 



statistic is the supremum of the absolute value of a Brownian Bridge. This immediately 
provides critical values without any bootstrap techniques. We impose a strong mixing 
assumption for the dependence structure. The proof relies on an invariance principle for 



multivariate sequential empirical processes from Biicher and Volgushev (2011). 



By using the copula-based expression for Spearman's rho from Schmid and Schmidt (2007) 



or 



Nelsen (2006), we get quite another contribution with our test, i.e. an extension of 



the copula constancy tests proposed by Busetti and Harvey (2011) and Kramer and van 



Kampen (2011). Since copula models are frequently used in financial econometrics (see 



e.g. Manner and Reznikova 2011 and Giacomini et al. 2009), such tests for structural 
change are important in this area. However, they are restricted to the case of testing 
for copula constancy in one particular quantile, e.g. the 0.95-quantile. This might be an 
important null hypothesis as well, but our test now allows for testing constancy of the 
whole copula by integrating over it. We therefore reject the null hypothesis of constant 
Spearman's rho (which is closely connected to the null hypothesis of an overall constant 
copula) if the integral over it fluctuates too much over time. The problem of testing 



constancy of the whole copula has recently be dealt with in the literature (see Ruppert 



2011, van Kampen and Wied, 2012 and Remillard 2010), but all these approaches need 



the computationally intensive bootstrap for approximating the limit distribution. With 
our approach, we need much less computational time for calculating the critical values. 
The paper is organized as follows: Section 2 presents our test statistic and the asymptotic 
null distribution, Section 3 considers local power, Section 4 presents Monte Carlo evidence 
about the test's behavior in small samples and Section 5 compares our new test with the 
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BPC test in terms of robustness by a simulation study and an empirical application. 
Finally, Section 6 concludes. All proofs are the appendix. 



2. Test statistic and its asymptotic null distribution 

In this section, we present the test statistic and the limit distribution of our test under the 
null. First, we introduce some notation: (Xi, . . . , X n ) are <i-dimensional random vectors 
on the probability space (Q, 21, P) with X,,- = (Xij, . . . , X d j). Regarding the dependence 
structure, we impose the following assumption: 

(Al) Xi, . . . , X„ are a- mixing with mixing coefficients atj satisfying 



OG 



^■2^7/(4+7) < ^ 



for some 7 G (0, 2) 



This dependence assumption is similar to the assumption made in Inoue (2001 ) and holds 



in most econometric models relevant in practice, e.g. for ARMA- and GARCH-process 



under mild additional conditions, see e.g. Carrasco and Chen (2002). 



The vectors Xy, j = 1, . . . , n, have joint distribution functions with 

i^'(x) = P(X 1)j <x u ..., x dij < x d ), x = On, . . . , x d ) G R d , 

and marginal distribution functions F it j(x) = P(Ajj < x) for x G 1R and i = 1, . . . ,d 
which are assumed to be continuous. 

According to Sklar's (1959) theorem, there exists a unique copula function Cj : [0, l] d — > 
[0, 1] of Xj with 

F*(x) = C j (F 1>j {x 1 ),...,F d j{x d )) 
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and 



Cj{u) = Fj(F{j(ui), F£(u d )), u = ( Ul , ...,u d )E [0, 1] £ 



where F 1 is the generalized inverse function, see e.g. Schmid and Schmidt (2007) 
In terms of the copula, Spearman's rho is defined as 



pj = h(d) ■ 2 



[2 d [ Cj(u)du- 1 j 



with h(d) = 2d d +J +1 ) , 



sec 



Schmid and Schmidt 



(2007) or 



Nelsen 



(2006). There are 



also other possibilities to define Spearman's rho in higher dimensions (see Schmid and 



Schmidt, 2007 and Quessy, 2009), but we focus on this expression for ease of exposition 



and because this measure performs well in terms of the asymptotic relative efficiency 



compared with other multivariate extensions, see Quessy (2009), p. 328. 
We are testing 



H : pj = p , j = 1, • • • , n vs. H x : 3j 6 {1, . . . , n - 1} : pj ^ p j+1 . 



Let 



1 

Fi. n {x) = - y j l{ Xi<j <x}, i = 1, • • • , d, x 6 R, 



U id = Fij(X itj ) and 



Uij-n ■= F Mn (Xij) = - ■ (rank of X itj in X ijU . . . , X i>n ) , i = 1, . . . , d, j = 1, . . . , n. 

71/ 
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Let R 3 (u) = l{u 1 , J <u 1 ,...,u d ,,<u d } and R^u) = ^ { u 1J;n < Ul ,...,u dJ;n <n d y The co P ula C is 
estimated by the empirical copula, defined as 



n d 



C n (u) = \ £ R 3 (n) = \ E II V^**)' U = -.«««) e [0, 1]' 

i=l j'=l i=l 



The estimator based on the first k observations is 



k d 



^( u ) = \ e ^( u ) = \ e n v^*,}' u = • • • > «*) g [°> ^ 



j=i »=i 



Note that we must use C^ijjn an d n °t to be able to apply the limit theorem from 



Biicher and Volgushev (2011) later on 



The estimator for the copula immediately yields an estimator for Spearman's rho: 



h = h(d) 



(V ^ ^ C k (u)du -lj= h{d) ■ J2 II^ 1 - Uiw) - 1 



We use the following test statistic W: 



W = D max 

Kk<n 



n 



{pk ~ Pn) 



--: D sup \P n (s) 

sS[0,l] 



with P n (s) = (p[ns] — Pn) and with an estimator D = -y=, where 



j=i i=i \ j=i i=i 



i - a 



i,y,n) 



+2 



n— 1 / \ I n—m d 

m 



x7n/ \ — W 
m=l \ jf=l *=1 



sEIId-fc* 



j=i »=i 



The kernel fc(-) is contained in the class /C2 of Andrews (1991) which guarantees positive 



semi-definiteness (we can choose e.g. the Bartlett-kernel); the bandwidth 7 n is chosen 
such that 7 n = 0(77,2). 
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For our central theorem, we need two additional assumptions: 
(Bl) (Xx, • • • , X n ) is strictly stationary. 

(B2) The copula C and marginal distribution functions F^j = Fi,i = 1, . . .,d are con- 
tinuous. 



Note that Assumption (Bl)| is in line with other tests for structural change, see e.g. Inoue 



(2001). 



Theorem 1. Under Hq and Assumptions (Al), (Bl), (B2) 



W sup \B(s) 

s6[0,l] 



where B(s) is a one- dimensional Brownian Bridge. 

Theorem [T] allows for constructing an asymptotic test. The main tool for the proof which 



can be found in Appendix A is an invariance principle from Biicher and Volgushev (2011 ) 
for the multivariate sequential rank order process 



[ns\ 



3=1 



\ns\ 



n 



[ns] 



which was introduced by Riischendorf (1976). This limit condition does not require any 



smoothness assumption on the derivatives of the copula, see Segers (2011) for a detailed 
discussion of this issue. 

There exists an interesting relationship between our test for constancy of Spearman's rho 



and the copula constancy tests proposed by Busetti and Harvey (2011) and Kramer and 



van Kampen (2011): One can show (see Appendix B) that our test is as a functional 



of the multivariate r-quantics on which these copula constancy tests base. But, in fact, 
while the other tests examine if the copula in a particular quantile is constant, we can 



S 



test for constancy of the whole copula by integrating over it. Note that the integral of 
the copula is a particular functional of it which has nice properties as it e.g. leads to 
a simple limit distribution which is free of nuisance parameters and therefore needs no 
bootstrap approximations. The limit distributions are more involved and/or difficult to 



Ruppert 


( 


2011 


)■ 


van Kampen and Wied 



(2012 


) and 


Remillard 


(2010 



3. Local power 

This section considers the local power of our test. Since the copula function of the 
random vectors under consideration changes with n, we now operate with a triangular 
array (X™,. . . , X"), but we suppress the index n for ease of exposition. Let C(u) be 
a copula and let C'(s, u) be another copula with an additional index parameter s. We 
consider local alternatives of the type 



C;(u) 



5 



)C(u) 
n J 



5 r' (3 
C -,u 

n \n 



By choosing, say, C (s, u) = [1 — g(s)] C (u)+ g(s)C" (u) for some copula C" (■) and some 
function g(-) bounded by 1 we obtain the sequence of correlations 



Pi 



1 - 



8 fj 

g - 

n \n 



8 fj 
Po + —/=g - ] pa- 

'n \n 



Choosing e.g. g(s) = l{ s >i/2} would lead to local alternatives in which the copula changes 
after the middle of the sample. A continuous function g would lead to continuously 
changing copulas against which our test has power as well. 

To deduce limit results for the sequence of local alternatives Q, we need some more 
assumptions: 



(CI) The analogue mixing condition |(A1) holds for the triangular array. 



(C2) The marginal distribution functions of (Xi, . . . ,X n ) do not depend on j and are 
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continuous. 



(C3) The joint copula for the random vectors (X 1; . . . , X n ) with lag I, 
Ck(u lf u a ) := P(X y < /••,/(»!) V,, ; < 



is specified to 

(7,y(ui,u 2 ) = I 1 
5 



Cj(ui,u 2 ) H C, 



$ 2 h ( 3 3 + 1 



<n \ \/n 



i i — , , u l> u 2 

n \n n 

C{xx x )Cf (—,u 2 ) +C(u 2 )C' (^m 
\ n J \n 



with a constant 5 G (0,1]. In this equation C/(-, •) is the (assumed to be continuous) 
joint copula of some sequence of stationary random vectors & with lag Z, C(-) is the 
(assumed to be continuous) copula of £j. Analogously, C t (-, •,-,■) is the copula of 
some sequence of stationary random functions rji(-) with lag I, C (-, ■) is the copula 
of ??«(•)• The term C'(r, t) is continuous in t for all r and there is (r, r', t) such that 
C"(r,*) ^ C'{r',t). 



The slightly cumbersome Assumption |(C3) is similar to Assumption B in Inoue (2001) 
and is required for applying the limit theorem for the sequential empirical process under 



local alternatives and mixing conditions from Inoue (2001). Passing each element of v to 
positive infinity yields equation ((!]). 
With these assumptions, we get 



Theorem 2. Under Assumptions (CI) and (C3) 



W sup 

se[o,i] 



B(s) + 5Dh(d)2° 



C* (t, u) dtdu — s 

[o,i] d 3o 3[o,i] d J>\ 



C* (t, u) dtdu 



where D is the probability limit of D under the null hypothesis. 
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With this theorem and Anderson's Lemma we can deduce that the asymptotic level is 
always larger than or equal to a, see Andrews (1997) or Rothe and Wied (2012). 



4. Finite sample behavior 



We investigate the test's finite sample behavior and compare it to the BPC test by 
simulating the empirical size under the null hypothesis and the empirical power under 
various alternatives. In order to make the simulation study not too lengthy, we always 
use the Bartlett kernel and bandwidth [log(n)] both for our test and for the BPC test 



(similar to Wied et al. 2012) and we restrict ourselves to the Student copula which is 



one of the most common copulas for modeling financial returns according to Cherubini 



et al. (2004), p. 181. For serial dependence, we assume a bivariate AL4(l)-process 



X t = e t + Qtt-x with 



^01 = 0.3 ^ 



V 



6o = 0.2 



In this situation the e t ,t G Z, are independent and identically distributed, following the 
ti-distribution with shape matrix 



S 



1 q 



, \q\ < 1- 



(2) 



In this case, we do not have finite fourth moments (even no finite first moment), which 
are required for the BPC test, and the Bravais-Pearson correlation is not even defined. 
However, also in such a situation, it might be interesting to know if there is a struc- 
tural change in the dependence structure. The null hypothesis is that X t has con- 
stant correlation of po = 0.4. Additionally we consider six alternatives Al to A6, in 
which the correlation jumps after the middle of the sample from p = 0.4 to p\ = 
0.6, 0.8, 0.2, 0.0, —0.2, —0.4, —0.6, respectively. For the simulations, we generate realiza- 
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tions e , ex,..., e n/2 and e n/2 +i, . . . , e n with qi = pi\J = 0, l|j Table |l| reports 
rejection frequencies at the significance level a = 0.05 based on 1000 repetitions and 
sample sizes n = 500, 1000, 2000. 

-Table Q] here - 

We see that the size of our test is kept and that the empirical power increases with the 
magnitude and n. It is slightly higher for increasing than for decreasing correlations, 
what is important in risk management. The BPC test is not applicable at all, because it 
cannot distinguish between null hypothesis and alternative. This is an expected behavior 
considering that the asymptotic variance of the empirical correlation coefficient is an 
unbounded function of the fourth moments of the population distribution and that, on the 
other hand, Spearman's correlation coefficient is invariant under monotonely increasing, 
componentwise transformations and hence little effected by heavy tails. 
We get a partly similar result when using the setup from above with the ts-distribution, 
constant correlation 0.4 and by adding one heavy outlier of size, say, (40,-100) to the 
sample at time c • 500, c = 0.05, 0.1, . . . , 1. If the outlier comes late, the BPC test almost 
always rejects the null hypothesis; if it comes early, the test almost never rejects it. This 
is an odd behavior and makes the test unsuitable for this outlier scenario. Our new test 
always keeps the size, see Figure [T] 

-Figure [T] here - 

Repeating the experiments from Table [T] with the ^-distribution and under serial inde- 
pendence yields a better behavior for the BPC test, but the test does not keep its size. 
Our test keeps its size and its power is in general higher than the BPC power for de- 
creasing correlations, for higher sample size it is even comparable or higher for increasing 
correlations, see Table [2] These results are important for financial applications because 

the existence of finite variance of financial returns is doubtful, see e.g. Grabchak and 

2 The choice of the qi is due to the fact that with this, the Pearson correlation would be equal to 
Po resp. p\ if it existed. Spearman's rho lies then closely to these values as numerical approximations 
suggest. 
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Samorodnitsky (2010). Interestingly, the power of our test becomes lower for higher 
serial dependence, as additional simulations show. 

-Table [2] here - 

In a setup for distributions with lower tails and no outliers, the power of our copula-based 
test is rather low compared to the BPC test, see Table [3] for exemplary results for the 
^-distribution and serial independence. The efficiency becomes lower when the distribu- 
tion comes "closer" to the Gaussian distribution where the usual empirical correlation 
coefficient is the maximum likelihood, i.e. the most efficient estimator of the correlation. 
Note however that there is compelling empirical evidence that financial returns are not 
Gaussian distributed. In general, decreasing correlations are detected better than increas- 
ing correlations. Under serial dependence, the power of our test is lower, as additional 
simulations show. 

-Table [3] here - 

The power of our new test becomes higher when testing for a correlation change in more 
than two dimensions and when there is an equal change in every component, see Table [4] 
for exemplary results for the trivariate ^-distribution under serial independence. In this 
case, a direct comparison with the BPC test is not possible because this test is restricted 
to two dimensions. 

-Table H] here - 



5. Robustness 



5.1. Simulation evidence The two major advantages of our test compared to the 



BPC test proposed by Wied et al. (2012) are its applicability without any moment con- 



ditions at all (as compared to the existence of fourth moments for the BPC test) and its 
appealing robustness properties. The latter shall be visualized by an instructive exam- 
ple. Both proposed fluctuation tests mainly derive their robustness properties from the 
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respective properties of the underlying correlation measure. The robustness properties of 
Spearman's rho, along with several other correlation estimators, are studied in detail in 



Croux and Dehon (2010). 



We sample a path (xt)t=i,...,5ao °f length n = 500 of the bivariate MA(1) process X 4 = 
et + Oet-i, where the et, t G Z, are i.i.d. with a centered bivariate Gaussian distribution 
and covariance matrix S. We choose 9 = diag(0.3, 0.2) and S as in equation ^ with q 
such that po — 0.4. 

We add one mild outlier to the sample by setting, say, x 2 88 to (20, —50). We denote the 
resulting contaminated example by (x^) t= i n , where w indicates weak contamination. 
Figure [2] visualizes the process 

k 

bk = D—=(r k -r n ), k = l,...,n, 
>n 



where f& denotes the Bravais-Pearson correlation coefficient based on Xi, ...,Xfc, and D 



is a deviation estimator that is described in Wied et al. (2012) and scales the process 



such that (&[ns])«e[o,i] converges to a Brownian Bridge. 

-Figure [2] here - 

-Figure [3] here - 

The BPC test statistic is then sup 1<fc< „ \b k \. The grey line in Figure [2] corresponds to the 
uncontaminated sample {x.t)t=i,...,m the black line to (x^)^^...^. The single outlier has 
a dramatic effect on the Bravais-Pearson test statistic and, in this example, causes the 
null hypothesis to be rejected at the significance level 0.05. The results are similar if the 
outlier is placed at a different position and if other distributions for the innovations are 
used. 

Alternatively we create a strongly contaminated sample (x%)t=i,...,n by randomly placing 
10 outliers in the second half of the sample. Each outlier is of the form (y t , —yt), where 
y t is drawn from the uniform distribution on [—1000, —100] U [100, 1000]. Both, fraction 
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and size of the outliers in (x|) t=1> ... )n are about 10 times as large as in (xf)t=i,..., n - Figure 
[3] depicts the process 

- A; 

%l) k = D—={p h - p n ), k = l,...,n, 

once being computed from the uncontaminated sample (grey line) and once from the heav- 
ily corrupted sample (Xt)t=i,...,n (black line). We witness a slight distortion of (V>fc)fc=i,->« 
as a result of the contamination, but the location of the maximizing point as well as the 
decision of the test are unaffected. Similar results are obtained for other realizations. 

5.2. Empirical relevance This subsection shows that the outlier scenario described in 
the previous subsection might indeed be relevant for a practitioner who analyzes struc- 
tural changes in the dependence of assets. This can be exemplarily seen in the time 
period around the Black Monday, 19th October 1987, i.e. in the time from January 1985 
to December 1989 using daily return data from Datastream (yielding n = 1262). Both 
the Dow Jones Industrial Average and the Nasdaq Composite daily returns are extremely 
negative at this day and the absolute values of these returns are much higher than the 
other ones at this time. One day after, the Dow Jones return is positive again, while 
the Nasdaq return remains negative - both on a high level compared to the means and 
standard deviations of all days. Table [5] shows the exact values. 

-Table [5] here - 

These outliers are reflected in the BPC test statistic, see Figure |4j part (a), for a visu- 
alization of the process (&fc)fc=i,...,i262 an d the peak around the Black Monday. On 19th 
October, the successively estimated correlations become very high, but due to the be- 
havior one day after, they fall down immediately. Both phenomena together lead to the 
peakj^] 

3 This is not exactly the same situation as in Table 1. However, by similar simulations with two 
outliers as in this application we can reproduce the peak of Figure 3 (a) as well. Thus, the figures give 
two different examples of a bizarre behavior of the test statistic which are both due to the construction 
of the Pearson correlation coefficient. 
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Applying the BPC test gives a test statistic of 1.447 (p- value of 0.030) such that the null 
hypothesis of constant correlation is rejected on the 5%-level, but the test statistic would 
be much lower without this peak and the null would not be rejected. 
Our Spearman test statistic is not affected by this peak - see Figure |i| part (b) for a 
visualization of the process (V , fc)fc=i,...,i262 - and the test statistic is equal to 0.886 (p-value 
of 0.412). Therefore, one should probably conclude that the dependence structure did 
not change seriously after the Black Monday. Similar results were obtained for other time 
periods around 19th October 1987. 

-Figure [4] here - 



6. Discussion 

We have proposed a new test for constancy of Spearman's rho which is much more robust 



against outliers than the BPC test previously suggested by Wied et al. (2012). 

Our test also allows for testing if the whole copula of multivariate random vectors is 

constant, therefore extending former suggested copula constancy tests. 

It is a task for further research to extend the dependence structure to functionals of iid- 

or even of mixing processes to enlarge the class of models in which our test can operate. 



A. Appendix section 

Proof of Theorem [7] 

Denote with Z°°(R fc ) the function space of all bounded functions from IR fc to M. 
Consider first P n (s): 

P n (s) = H (h( d ) ■ (V ^ C [ns] (u)du -lj- h{d) • (V J ^ C n (u)du - 1^ 

r 1 p ( 1 M 1 n \ 



h(d) ■ 2 d / A n (s,u)du 

'[0,l] d 
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with 



A n (s,u) 



\ns\ 



n 



R E4-(«)-„E^) 



(3) 



With an invariance principle for 



^f>>)-C(u)), 

i=l 



(4) 



which is presented in Theorem 2.1 in Inoue ( 2001[ ) (setting S from Inoue, 2001 to and 



defining x ni from Inoue, 2001 as ,j))i<i<d), Condition 3.1 in Biicher and Volgushev 



(2011) is satisfied with G* n (s,u) := ^J2i=i R A n ) and c *( s i u ) sC ( u ) which lies the 



space defined in Biicher and Volgushev| ( 2011[ ). So we obtain with Corollary 3. 3. a in 



Biicher and Volgushev (2011) 



A n (-,-) ^ d Aq(;-) in I 



with 



A (s,u) = V (s,u) - sV (l,u). 

Vo(s,u) is a P-almost surely continuous, centered Gaussian process with covariance func- 
tion 



Kq^sx, ui), (s 2 , u 2 )) := Cov(Vb(si, ui), V (s 2 , u 2 )) = (si A s 2 )if'(ui, u 2 ), 



17 



where 



K'{n u u 2 ) = C( Ul A u a ) - C(u 1 )C(u 2 ) 

oo 

+ XI ( E ( 1 {X lj i<F 1 - 1 K) i ...;X dil 



<F d - 1 (^)} 1 {X ljm <F 1 - 1 («i);...;X (ilTO <F7 1 («^)} 



m=2 



E ( 1 {Xi,i<F- 1 K);...;X d , 1 <F d - 1 (nf)}) ' E ( 1 {X 1 , m <F- 1 ( M i);...;X d , m <F- 1 ( U J 



)}< 



m=2 

-E(l 



{Xi,!^^- 1 ^);...;^,!^^- 1 ^)}) ' E ( 1 {X 1 , m <F- 1 (^);...;X d , m <F- 1 («f)}, 



This covariance function is the limit of the covariance function of 



[ns] 

V n (s,u) = — V(^(u)-C(u 

7 = 1 



)) 



i.e. 



#o((si,ui),(s 2 ,u 2 )) = lim Cov(K(si,Ui),K(s 2 ,u 2 )), 



see Inoue (2001). Now, with the Continuous Mapping Theorem, 



where 



P ( s ) = h{d)2 d I A {s,u)du 
'[o,i] d 



is a P-almost surely continuous, centered Gaussian process. With Fubini's theorem, the 
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covariance function is 



Cov(P ( Sl ),P ( S2 )) 

= h{d) 2 2 2d [ [ Cov(y (si,ui) -s 1 V (l,u 1 ),V (s 2 ,u 2 ) - s 2 V (l,u 2 ))du 1 du 2 

J[0,l] d J[0,l] d 

= h(d) 2 2 2d ( Sl A 

s 2 ~~ s l s 2 ~~ s l s 2 + s 1 s 2 ) / / Cov(Vo(l,Ui),H(l,u 2 ))tiuitZu2 

J[0,l] d J[0,l] d 

= (si A s 2 - sis 2 )D' 



with 



D' = h(d) 2 2 2d [ [ K'(u 1 ,u 2 )du 1 du 2 

J[o,i] d J[o,i] d 

= h(d) 2 2 2d J E ( - tfy 



+2 



\i=i 

E E ( II( 1 - U *J)0- - U ^rn)\ ~ [ E ( 11(1 - C/ 5 



m=l \i=l 



.i=l 



= /l(d)2 



2d 



m=l \i=l 



i=l 



This holds because, again with Fubini, 



C(ui A u 2 )dui<iu2 

= /o i] d I[o ir E ( 1 ( x ^ F r 1 wy>----' x ^ F d 1 ^ 

= E U[o i] d I[o i] d 1 l x ^ F rH»{);---*dj<FzH*^^ 

= E ( / , / , 1 {^ij<«ii"-i^<«?} 1 {^ij<«5i"-i^<^} dUldu2 ) 
\./[o,i] d ^[o,i] d / 

/ d d 



,j=i 



i=l 



The other summands of K'(u 1: u 2 ) are integrated analoguesly. 
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We get a consistent estimator for D' from de Jong and Davidson (2000) 



n d 



n d 



& = Kdf^ a y: rid - u itj f - 1 1 x; na - ^ 



+2 



j=i i=i 

n— 1 / \ / n— m 

m 



i=l i=i 



n (i 



E * (7) E s IK 1 - *w - w - U E^ 1 - ^ 



n 

i=i »=i 



j=i i=i 



with a kernel that is contained in the class /C2 of Andrews (1991) which guarantees 



positive semi-defmiteness of D'. Next, we show that D' — D' —¥ p 0. By the invariance 
principle Q (this time applied on the components of Xj) we get a Gliwenko-Cantelli-like 
theorem (in probability) with rate for the marginal empirical distribution functions, 
that means, 



B n : = max sup \ Uij ]n — U^j] = Op In 2 



=l,...,d 



Since < Uij, Uij- n < 1, we obtain 



l[(l-U i}j . n )(l-U, 



i,j+m;n ) 



i=l 



t=l 



< 2dB„ 



Thus we get 



n-l , x 

D' - D| < C X A; B n = P ( 7n rH 



o P (l) 



Andrews 



(1991), p. 852. Therefore D' is a 



as 7 n = 0(712), compare the argument in 
consistent estimator of D' . 

The theorem follows then with the Continuous Mapping Theorem, because the process 



D> 



is a P-almost surely continuous, centered Gaussian process with the same covariance 
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function as the Brownian Bridge, i.e. 

Cov(P *(si), P * (s 2 )) = S 1 AS 2 ~ S X S2- 

Since a Gaussian process is uniquely determined by the first two moments, the limit pro- 
cess is in fact a Brownian Bridge. ■ 



Proof of Theorem [2| 

The proof is basically similar to the proof of Theorem [l] We consider A n (s,u) from 
(|3| and use a minor modification of Corollary 3. 3. a in |Biicher and Volgushev| (2011). 
Condition 3.1 is modified in that sense that we formulate the condition for triangular 



arrays. The invariance principle from Theorem 2.1 in Inoue (2001) (defining again x. 



from Inoue, 2001 as (Fi(Xi,j))i<i<d with the stationarity of the marginal distribution 
functions) yields that 



[ns\ 

' it , 



converges in distribution against the process Vq(s,u) such that 



(5) 



\TIS I [?"ISJ |^TiSj 

)= - C{u)) = -L ^(i?,(u) - Cj(u)) + 4= £(G,(u) " C ^ ( 6 ) 

v 7=1 



i=l 



|7isj [fis 

J>i(u) - C,(u)) + -J2( C '( J -^)~ °W ) ( 7 ) 



converges in distribution against the process Vq(s, u) + 5 f* C (t, u) dt — SsC(u). This pro- 

Biicher and Volgushev| (p01l\ with G* n (s, u) = \ R j( u ) 



cess fulfills Condition 3.1 in 



and C*(s, u) = sC(u) which lies the space defined in Biicher and Volgushev (2011). 
Thus the proof of Corollary 3. 3. a goes through in the same way as in the stationarity 
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case. Thus, A n (s,u) converges against 



A (s,u) + 5 



C (t,u)dt-s C (t, u) dt 
Uo Jo 



In addition, the probability limit of D under the sequence of local alternatives is the 
quantity D from the proof of Theorem [lj i.e. the probability limit of D under the null 
hypothesis. This holds because 

lim Cov(y n (si,Ui),K(s 2 ,u 2 )) 



is the same under the null hypothesis as well as under the sequence of local alternatives. 
Thus, the theorem is proved. ■ 



B. Connection to copula constancy tests 



Let F i ^(ui) denote the empirical quantile function. Write 



AJs,u) 



1 / 1 L nS > 1 n 



j'=i 



n \ [nsl ^— ' 

\ L J 5=1 

I '"^ I L ' / 1 " 



J'=l 

[ns] 



/n ^ 



j'=i 

[ns] 



C n {u) - Rj(u] 



k=l \ j=l 
1 n 



?? 



m=l 



-= £ - I(X y < F^( Ul ), ...,X dj < F^(u d )) 

.7=1 
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where the last step uses that l/"^X]m=i Rm(u) = C n {u) and 



Rj{u) 



I{U\j,n < "1, • • • , Udj,n < U d ) 

nKn(Ulj,n) < *g(t*l), F^(U dj , n ) < F^{u d )) 

I(X lj <F^(u 1 ),...,X cij <F^(u d )). 



Note that C n {u) — I{X\j < F l7 ^(ui), ■ ■ ■ -,X d j < F d ^{u d )) are the bivariate r-quantics in 



Busetti and Harvey (2011) if d = 2. Hence, the test can be written as an functional form 



of the same quantities. 



References 

Amaral, L., V. Plerou, P. Gopikrishan, M. Meyer, and H. Stanley (2000): 
"The distribution of returns of stock prices," International Journal of Theoretical and 
Applied Finance, 3(3), 365-369. 

Andrews, D. (1991): "Heteroskedasticity and autocorrelation consistent covariance 
matrix estimation," Econometrica, 59(3), 817-858. 

- (1997): "A conditional Kolmogorov test," Econometrica, 65, 1097-1128. 

Aue, A., S. Hormann, L. Horvath, and M. Reimherr (2009): "Break detection in 
the covariance structure of multivariate time series models," The Annals of Statistics, 
37(6B), 4046-4087. 

Bucher, A. AND S. Volgushev (2011): "Empirical and sequential empirical copula 
processes under serial dependence," preprint. 

Busetti, F. and A. Harvey (2011): "When is a copula constant? A test for changing 
relationships," Journal of Financial Econometrics, 9(1), 106-131. 

CARRASCO, M. AND X. CHEN (2002): "Mixing and moment properties of various 
GARCH and stochastic volatility models," Econometric Theory, 18, 17-39. 

23 



Cherubim, U., E. Luciano, and W. Vecchiato (2004): Copula methods in finance, 
Wiley. 

CROUX, C. AND C. DEHON (2010): "Influence functions of the Spearman and Kendall 
Correlation measures," Statistical Methods and Applications, 19(4), 497-515. 

de Jong, R. and J. Davidson (2000): "Consistency of kernel estimators of het- 
eroscedastic and autocorrelated covariance matrices," Econometrica, 68(2), 407-424. 

Embrechts, P., A. McNeil, and D. Straumann (2002): "Correlation and depen- 
dence in risk management: Properties and pitfalls," in Risk management: Value at Risk 
and beyond, ed. by M. Dempster, Cambridge University Press, Cambridge, 176-223. 

Gaissler, S. and F. Schmid (2010): "On testing equality of pairwise rank correlations 
in a multivariate random vector," Journal of Multivariate Analysis, 101, 2598-2615. 

Galeano, P. (2007): "The use of cumulative sums for detection of changepoints in 
the rate parameter of a Poisson Process," Computational Statistics and Data Analysis, 
51(12), 6151-6165. 

Galeano, P. and D. Pena (2007): "Covariance changes detection in multivariate time 
series," Journal of Statistical Planning and Inference, 137(1), 194-211. 

Giacomini, E., W. Hardle, and V. Spokoiny (2009): "Inhomogeneous Dependence 
Modeling with Time- Varying Copulae," Journal of Business and Economic Statistics, 
27(2), 224-234. 

GRABCHAK, M. AND G. SAMORODNITSKY (2010): "Do financial returns have finite 
or infinite variance? A paradox and an explanation," Quantitative Finance, 10(8), 
883-893. 

Inoue, A. (2001): "Testing for distributional change in time series," Econometric The- 
ory, 17, 156-187. 



24 



KRISHAN, C, R. PETKOVA, AND P. RlTCHKEN (2009): "Correlation risk," Journal of 
Empirical Finance, 16, 353-367. 

Kramer, W. (2002): "Statistische Besonderheiten von Finanzzeitreihen," Jahrbiicher 
fur Nationalokonomie und Statistik, 222, 210-229. 

Kramer, W. and M. van Kampen (2011): "A simple nonparametric test for structural 
change in joint tail probabilities," Economics Letters, 110, 245-247. 

Lee, S., J. Ha, O. Na, and S. Na (2003): "The CUSUM test for parameter change 
in time series models," Scandinavian Journal of Statistics, 30(4), 781-796. 

LONGIN, F. AND B. SOLNIK (1995): "Is the correlation in international equity returns 
constant: 1960-1990?" International Money and Finance, 14(1), 3-26. 

MANNER, H. AND O. REZNIKOVA (2011): "A survey on time-varying copulas: Specifi- 
cation, simulations and application," Econometric Reviews, 31(6), 654-687. 

Nelsen, R. (2006): An introduction to copulas, Springer, 2 ed. 

Ploberger, W., W. Kramer, and K. Kontrus (1989): "A new test for structural 
stability in the linear regression model," Journal of Econometrics, 40, 307-318. 

Quessy, J.-F. (2009): "Theoretical efficiency comparisons of independence tests based 
on multivariate versions of Spearman's rho," Metrika, 70, 315-338. 

Remillard, B. (2010): "Goodness-of-fit tests for copulas of multivariate time series," 
preprint. 

Rothe, C. and D. Wied (2012): "Misspecification testing in a class of conditional 
distributional models," preprint. 

RUPPERT, M. (2011): "Consistent testing for a constant copula under strong mixing 
based on the tapered block multiplier technique," preprint. 



25 



Ruschendorf, L. (1976): "Asymptotic distributions of multivariate rank order statis- 
tics," The Annals of Statistics, 4(5), 912-923. 

Schmid, F. AND R. Schmidt (2007): "Multivariate extensions of Spearman's rho and 
related statistics," Statistics and Probability letters, 77, 407-416. 

Segers, J. (2011): "Asymptotics of empirical copula processes under nonrestrictive 
smoothness assumptions," Bernoulli, forthcoming. 

Sklar, A. (1959): "Fonctions de repartition a n dimensions et leurs marges," Publica- 
tions de I'Institut de Statistique de L'Universite de Paris, 8, 229-231. 

VAN KAMPEN, M. AND D. WlED (2012): "A nonparametric constancy test for copulas 
under mixing conditions," preprint. 

WlED, D., W. Kramer, and H. Dehling (2012): "Testing for a change in correlation 
at an unknown point in time using an extended functional delta method," Econometric 
Theory, 58(3), 570-589. 



26 




0.2 0.4 0.6 0.8 1.0 

Timing outlier point 



Figure 1: Empirical rejection frequencies in an outlier scenario 



Table 1: Empirical power in different settings for the ti-distribution, po = 0.4, results of the 
BPC test in brackets 



n 








Values of p\ 










0.4 


0.6 


0.8 


0.2 





-0.2 


-0.4 


-0.6 


500 


0.044 


0.067 


0.165 


0.071 


0.138 


0.282 


0.506 


0.740 




(0.489) 


(0.514) 


(0.560) 


(0.480) 


(0.481) 


(0.502) 


(0.523) 


(0.540) 


1000 


0.048 


0.102 


0.282 


0.094 


0.256 


0.530 


0.811 


0.969 




(0.479) 


(0.512) 


(0.552) 


(0.474) 


(0.496) 


(0.497) 


(0.514) 


(0.544) 


2000 


0.051 


0.160 


0.486 


0.155 


0.466 


0.834 


0.982 


1.000 




(0.487) 


(0.495) 


(0.557) 


(0.486) 


(0.467) 


(0.500) 


(0.511) 


(0.538) 
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0~ 100 200 300 400 500 



Figure 2: Process (bk)k=i,...,500 computed from weakly contaminated (black) and uncontami- 
nated (grey) sample 




-0.5 

6 100 200 300 400 500 



Figure 3: Process (V ; fc)fc=i,...,500 computed from strongly contaminated (black) and uncontami- 
nated (grey) sample 
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Table 2: Empirical power in different settings for the i2-distribution, po = 0.4, results of the 
BPC test in brackets 



n 








Values of p 1 










0.4 


0.6 


0.8 


0.2 





-0.2 


-0.4 


-0.6 


500 


0.037 


0.073 


0.224 


0.083 


0.216 


0.449 


0.736 


0.951 




(0.132) 


(0.243) 


(0.486) 


(0.128) 


(0.242) 


(0.387) 


(0.513) 


(0.633) 


1000 


0.050 


0.122 


0.426 


0.130 


0.394 


0.783 


0.974 


0.999 




(0.134) 


(0.220) 


(0.508) 


(0.146) 


(0.259) 


(0.444) 


(0.567) 


(0.720) 


2000 


0.049 


0.229 


0.697 


0.203 


0.710 


0.983 


1.000 


1.000 




(0.115) 


(0.249) 


(0.580) 


(0.148) 


(0.301) 


(0.478) 


(0.664) 


(0.745) 



Table 3: Empirical power in different settings for the t3-distribution, po = 0.4, results of the 
BPC test in brackets 



n Values of pi 





0.4 


0.6 


0.8 


0.2 





-0.2 


-0.4 


-0.6 


500 


0.045 


0.085 


0.212 


0.084 


0.248 


0.506 


0.778 


0.951 




(0.057) 


(0.288) 


(0.739) 


(0.110) 


(0.400) 


(0.695) 


(0.853) 


(0.904) 


1000 


0.043 


0.124 


0.385 


0.138 


0.441 


0.782 


0.973 


1.000 




(0.044) 


(0.349) 


(0.856) 


(0.185) 


(0.559) 


(0.843) 


(0.914) 


(0.951) 


2000 


0.051 


0.239 


0.700 


0.212 


0.715 


0.986 


0.999 


1.000 




(0.040) 


(0.487) 


(0.895) 


(0.294) 


(0.766) 


(0.900) 


(0.944) 


(0.962) 



Table 4: Empirical power in different settings for the trivariate t3-distribution, po = 0.4 (for 
pi = —0.6, the shape matrix would not be positively definite) 



n 






Values of 


Pi 








0.4 


0.6 


0.8 0.2 





-0.2 


-0.4 


500 


0.040 


0.154 


0.545 0.206 


0.687 


0.988 


1.000 


1000 


0.047 


0.310 


0.850 0.381 


0.950 


1.000 


1.000 


2000 


0.042 


0.574 


0.996 0.637 


0.999 


1.000 


1.000 



Table 5: Dow Jones and Nasdaq Returns around the Black Monday and empirical moments of 
the whole time span 



Day 


16.10. 


19.10. 


20.10. 


21.10. 


22.10. 


p 


a 


Dow Jones 


-0.047 


-0.256 


0.057 


0.097 


-0.039 


0.000 


0.013 


Nasdaq 


-0.039 


-0.120 


-0.094 


0.071 


-0.046 


0.000 


0.009 
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